home *** CD-ROM | disk | FTP | other *** search
/ Languguage OS 2 / Languguage OS II Version 10-94 (Knowledge Media)(1994).ISO / gnu / libg_261.zip / libg_261 / libg++ / src / pow.cc < prev    next >
C/C++ Source or Header  |  1992-01-17  |  1KB  |  71 lines
  1. /* 
  2. Copyright (C) 1990 Free Software Foundation
  3.     written by Doug Lea (dl@rocky.oswego.edu)
  4.  
  5. This file is part of the GNU C++ Library.  This library is free
  6. software; you can redistribute it and/or modify it under the terms of
  7. the GNU Library General Public License as published by the Free
  8. Software Foundation; either version 2 of the License, or (at your
  9. option) any later version.  This library is distributed in the hope
  10. that it will be useful, but WITHOUT ANY WARRANTY; without even the
  11. implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
  12. PURPOSE.  See the GNU Library General Public License for more details.
  13. You should have received a copy of the GNU Library General Public
  14. License along with this library; if not, write to the Free Software
  15. Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.
  16. */
  17.  
  18. #ifdef __GNUG__
  19. #pragma implementation
  20. #endif
  21. #include <builtin.h>
  22. #include <math.h>
  23.  
  24. double pow(double x, long p)
  25. {
  26.   if (p == 0)
  27.     return 1.0;
  28.   else if (x == 0.0)
  29.     return 0.0;
  30.   else
  31.   {
  32.     if (p < 0)
  33.     {
  34.       p = -p;
  35.       x = 1.0 / x;
  36.     }
  37.  
  38.     double r = 1.0;
  39.     for(;;)
  40.     {
  41.       if (p & 1)
  42.         r *= x;
  43.       if ((p >>= 1) == 0)
  44.         return r;
  45.       else
  46.         x *= x;
  47.     }
  48.   }
  49. }
  50.  
  51. long  pow(long  x, long p)
  52. {
  53.   if (p == 0)
  54.     return 1;
  55.   else if (p < 0 || x == 0)
  56.     return 0;
  57.   else
  58.   {
  59.     long r = 1;
  60.     for(;;)
  61.     {
  62.       if (p & 1)
  63.         r *= x;
  64.       if ((p >>= 1) == 0)
  65.         return r;
  66.       else
  67.         x *= x;
  68.     }
  69.   }
  70. }
  71.